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ANALYSIS OF HIGH-ASPECT-RATIO 
JET-FLAP WINGS OF ARBITRARY GEOMETRY 

By Peter B. S. Lissaman 
Northrop Corporate Laboratories 


SUMMARY 

Many techniques of amplifying the lift of a wing involve the blowing of 
a thin, high energy sheet from the wing trailing edge (jet flap). Currently, 
no methods exist to rapidly compute the performance of an arbitrary 
jet- flapped wing. 

This paper develops an analytical technique based on the method of Maskell 
and Spence in which the well-known lifting-line approach is coupled with an 
auxiliary equation providing the extra function needed in jet-flap theory. 

The method of Maskell and Spence was developed for a special class of wings 
of elliptical planform and loading. The present method generalizes this to 
handle straight, uncambered wings of arbitrary planform, twist, and blowing 
(including unsymmetrical cases). An analytical procedure is developed for 
continuous variations in the above geometric data with special functions to 
exactly treat discontinuities in any of the geometric and blowing data. A 
rational theory for the effect of finite wing thickness is introduced as well 
as simplified concepts of effective aspect ratio for rapidhand estimation 
of performance. 

The procedure has been programed and prints out all relevant data. 

Each wing case takes less than six seconds on the CDC 6600. 



The results converge properly to all known limit cases, and correlation 
with other known theoretical solutions for jet-flap wings is good. It is found 
that, although the method is a lifting-line technique, good accuracy was 
obtained for aspect ratios as low as 4. Correlation with experimental data 
was as good as other techniques, but it is pointed out that few properly 
controlled and corrected experimental tests are known. 

It is believed that the accuracy of this method is quite adequate for 
jet-flap wing design and that the speed and flexibility of the program make 
it very suitable for preliminary design. 

INTRODUCTION 

Powered lift is becoming an increasingly important technique for STOL 
vehicles. In general, this consists of methods in which the energy of a 
certain mass of air is increased, and this mass then forced over aerodynamic 
surfaces to amplify the lifting effects. Typical examples of powered lift 
using a jet sheet with lifting-surface interaction are: jet-flap wing, 

externally blown flap, and augmentor wing. The theory developed in this 
analysis is applicable to all these types of lift amplification, since augmentor 
wings and externally blown flaps involve a similar analysis to jet-flap wings 
and may be approximated by thin sheet jet- flap theory. 

The fluid dynamics of multi -energy flows is generally a poorly understood 
field. Although it is simple to write the basic equations satisfied at appro- 
priate points, a full solution for a specified geometry is more difficult, even 
neglecting the effect of viscosity. In essence, this problem occurs because 
the unknown boundary between the two flows of different energy cannot be specified 
in advance. Since this boundary moves as the flow state is altered, the solutions 
are intrinsically nonlinear because of the nonlinear boundary conditions. 

For the problem considered, a high energy jet sheet issuing from the 
trailing edge of the wing is assumed. The standard assumptions of incompressible 
inviscid flow, a thin jet, and linearized boundary conditions are made. This 
model is described in detail in reference 1. 
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SYMBOLS 


a 


A 

A. 


c 

c' 

c* 



C 

C 


£ 


intermediate complex mapping parameter 
h 2 

aspect ratio, £_ 

effective aspect ratio (see eq. (59)) 

coefficients of $ c sine expansion 

wing span, meters (ft) 

local wing chord, meters (ft) 

chord of airfoil in auxiliary plane (app. B) 

normalized chord, 

c av 

s 

mean chord, , meters (ft) 

b 

effective chord, (see eq. (19)) 

1 / rrr 

effective chord at station, oj = y- 

momentum coefficient in auxiliary plane (app. B) 

momentum coefficient in auxiliary plane at trailing edge 
(see eq. (B12)) 

chord of airfoil in physical plane, meters (ft) (see app. B) 
lift coefficient at a station 
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local Tift coefficient due to jet Inomentum 

two-dimensional normal pressure lift coefficient (see eq. (6)) 

two-dimensional lift coefficient finite thickness airfoil 
two-dimensional lift coefficient zero thickness airfoil 
two-dimensional lift coefficient gradient with angle of attack 

two-dimensional lift coefficient gradient with jet-flap deflection angle 

local induced-drag coefficient 

local induced-drag coefficient due to jet momentum 

local induced-drag coefficient due to pressure on airfoil 

wing- induced drag coefficient 

local jet-momentum coefficient, ^ 

\ P V 2 c 

wing jet-momentum coefficient 
wing lift coefficient 

wing lift coefficient gradient with angle of attack 
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wing lift coefficient gradient with jet angle 


wing pressure lift coefficient 

wing pressure lift coefficient gradient with angle of attack 

wing pressure lift coefficient gradient with jet angle 

local normal -force coefficient in airfoil 
local nose-thrust coefficient on airfoil 
induced drag of wing 

elements of matrix determining a change 
induced drag efficiency 

function used to define jump operator (see eq. (37)) 
ratio of blown to unblown lift slope of airfoil (see eq. (59)) 
mapping derivative from physical to auxiliary plane 
elements of line vector denoting error in nose thrust 

components of force in Trefftz plane (app. A) 


Heaviside operator 



L 

m,n,i ,j 
N 




q 


elements of matrix defining discontinuous potential gradients 

lift of wing, grams (lbs) 
indices of summation 
leading edge of airfoil (app. B) 
nose source 

gradient of nose source with angle of attack 

gradient of nose source with jet angle 

complex conjugate velocity in auxiliary plane 
free-stream speed in auxiliary plane, meters/second (ft/sec) 


Q 


r 


R 

s 

S 


complex conjugate velocity in physical plane 
free-stream speed in physical plane, meters/second (ft/sec) 

even integer defining upper limit of potential function series, 
(see eq. (22)) 

radius of curvature of wake in physical plane 
arc length in auxiliary plane (app. B) 
arc length in physical plane (app. B) 
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2 2 

S wing area, meters (ft ) 

w 

t thickness of airfoil in auxiliary plane (app. B) 

T thickness of airfoil in physical plane (app. B) 

T' trailing edge of airfoil (app. B) 

V main-stream velocity, meters/second (ft/sec) 

w(z) normalized complex perturbation potential in Trefftz plane 

W complex potential function 

x spanwise coordinate in Trefftz plane 

X spanwise coordinate normalized by semi span (2x/b) 

y vertical coordinate in Trefftz plane 

Y vertical coordinate normalized by semispan (2y/b) 

z complex coordinate in Trefftz plane (app. A) 

z' complex coordinate in auxiliary plane (app. B) 

Z complex coordinate in physical plane (app. B) 

a local wing angle of attack, radians 

a mean angle of wake in Trefftz plane, radians 

a e effective angle of attack, radians (see eq. (20)) 
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a 


e 


k 




a 


2 


a 


... , _ |(*jy . 

effective angle of attack at station, “ = ~ , radians 
local induced angle of attack, radians 

local effective angle of attack of jet-flap airfoil, radians, 
(see eq. (13)) 

local effective angle of attack of plain airfoil, radians 
local angle of downwash in Trefftz plane, radians 


^kn 

3*(w) 
m v ' 


local angle of jet wake in Trefftz plane, radians 
elements of matrix defining potential gradient 

elements of line vector defining potential gradient at uj 


e**( w 



y(x) 

Y a 


r 



elements of line vector used in discontinuous analysis 

local jet-moment per unit span at station x 
circulation about airfoil in auxiliary plane (app. B) 

circulation about airfoil and wake in physical plane, (app. B) 
circulation about airfoil alone in physical plane (app. B) 
Kronecker delta 


Ap 


pressure perturbation in Trefftz plane (app. A) 
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e 

n 


0 

e 

e 

x. 

u 

p 


a 


$ 


$ 


c 


<$>* 


V*y 

x o 


a) 



constant used in iteration technique (see eq. (47)) 

ratio of lift slope of thick unblown airfoil to zero thickness 
unblown airfoil (see eq. (54)) 

local jet angle relative to chordline, radians 

slope of wake in auxiliary plane (app. B) 

slope of wake in physical plane (app. B) 

magnitude of i^* 1 downwash discontinuity (see eq. (35)) 

intermediate complex mapping parameter 

3 3 

density of main-stream flow, grams/meter (slugs/ft ) 

function relating Trefftz plane downwash to wing plane downwash 

normalized potential function 

continuous normalized potential function 

discontinuous normalized potential function 

gradient of potential function, w.r.t. x,y 

distance from leading edge of airfoil 

angular coordinate of wing station (cos - ^ x) 

x L 

angular coordinate of i in discontinuity of wing geometry 
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- - PREVIOUS WORK 

The linearized problem for the two-dimensional jet flap both in and out 
of ground effect has been treated in references 2 to 5. Excellent correlation 
of the theory with experiment has been found. 

The three-dimensional problem has been handled by a number of workers 
using different approaches to handle the kinematic and dynamic boundary condi- 
tions in the wake. A semi -empirical method for computing wing performance is 
given by Kuchemann (ref. 6). Solutions using the techniques of matched 
asymptotic expansions have been published by Kerney (ref. 7) and Tokuda (ref. 8). 
All of these methods involve approximations about the wake location. 

The model by Maskell and Spence (ref. 1) makes an attempt to handle the 
dynamic condition on the wake by using two-dimensional solutions which properly 
match the wake boundary conditions. However, this treats only the high-aspect- 
ratio elliptically loaded case. Perhaps the most advanced model is that given 
by Lopez and Shen (ref. 9). Although this model (like the others) is linearized, 
it utilizes lifting-surface theory, finding the wake deformation by iteration 
on the proper dynamic and kinematic boundary conditions. 

The model proposed here lies somewhere between the method of Lopez and 
Shen (ref. 9) and that of Maskell and Spence (ref. 1) in its assumptions and 
complexity. It is essentially an extension of the blown lifting-line concept 
for a high-aspect-ratio wing of arbitrary planform and blowing. In the current 
model, methods are developed to take into account discontinuities in the wing 
geometry in chord, angle of attack, jet angle, or blowing coefficient. 

ANALYSIS 

General 

The formal lifting-surface solution consists of matching kinematic boundary 
conditions on the wing as well as kinematic and dynamic boundary conditions on 
the loaded jet wake as shown by Lopez and Shen (ref. 9). An extension of the 
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Prandtl lifting-line approach as employed by Maskell and Spence (ref. 1) 
is adopted. 

The basic idea (which may be formalized in terms of matched asymptotic 
expansions) is that the forces on the wing may be determined in two places. 

In the Trefftz plane, the forces may be uniquely determined in terms of wake 
parameters (essentially downwash and jet coefficient). On the wing itself 
(wing plane), the forces may be determined in terms of local section properties. 
Here, the main parameters are the blown airfoil characteristics and the induced 
downwash angle at the wing. The Trefftz plane and wing plane forces are then 
matched to provide sufficient equations for solution. 

In the case of an unblown wing, the lift in the two planes is matched to 
provide a relationship between the potential <f> and the Trefftz plane down- 
wash <j> generating the well-known integro-differential equation of lifting- 
line theory (Prandtl's equation). An important assumption is that the 
induced angle at the wing is taken to be half that in the Trefftz plane. 

In the case of a blown wing, a major portion of the bound vorticity 
(that directly due to the jet sheet) is located downstream of the trailing 
edge. Thus, it is not valid to assume that the trailing vorticity emanates 
from a lifting-line situated within the wing chord. 

The analysis involves only two space variables, with the coordinate system 
centered on the wing or wake centerline. The independent variable x is used 
for the spanwise dimension, positive towards the right wingtip where the wing 
is viewed from the rear, with y the vertical coordinate positive in the 
upwards (lift) direction. (See fig. 1.) 

The approach of reference 1 is adopted and introduces a new constant a 
(variable across the span) by which the Trefftz plane downwash is scaled. 

Thus, the downwash at the wing is assumed constant across each chordwise 
station and equal to a(x) a oo (x) where a oo (x) is the downwash angle in 
the Trefftz plane. Hence, the induced angle a..(x), as shown in figure 2, is 
given by 

a.,- (x) = o(x) a^x) (1) 
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This implies that a further unknown is introduced into the wing plane 
solution, so that an additional matching equation is required. It seems natural 
then to match both the lift and drag equations, or equivalently, the normal and 
chordwise force equations providing the additional equation for a . Reference 1 
assumes that a is constant across the span, as implied by the elliptical loading. 
In the present solution, a is permitted to vary with x and determined at each 
spanwise station. 


Trefftz Plane Solution 

By considering the flow in the Trefftz plane and by resolving parallel and 
normal to the chord, a pair of equations is obtained for the normal pressure 
coefficient and the thrust coefficient on the wing. It should be noted that 
the total force on the wing constitutes the sum of the pressure force and the 
momentum force due to the jet. Since the latter can be directly computed from 
the kinematics, it is more convenient to work with the pressure forces alone. 

It is noted that all coefficients are normalized by division of the dimensional 

2 2 

force by 1/2 P V c for the two-dimensional quantities and by 1/2 p V S w 

for the integrated wing quantities, where P is the ambient density , V the 
free-stream velocity, c the local chord, and S w the wing area. Appendix A 
shows these force coefficients to be given by 

C n ’ *4* ' C j ( 5 +y + “ + e ) 

• ••,)'} 

Wing Plane Solution 

This solution is required to relate the local wing angle of attack to its 
bound vorticity distribution. This distribution must be such that the downwash 
it causes when added to the trailing system flow produces the downwash necessary 
to match the physical kinematic boundary conditions on the wing (fig. 2). 


( 2 ) 

(3) 
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Consistent with lifting-line theory, it is assumed that the local downwash 
may be computed by considering the downwash of a two-dimensional system as shown 
in reference 1. This solution must match the kinematic and dynamic boundary 
conditions appropriate to a jet-flapped airfoil. 

There are two approximate approaches to the problem of defining the 
boundary vorticity. The first is to assume that the vorticity is composed of 
two parts - the one being that of a two-dimensional jet-flap airfoil at 
a l = a " °oo with a jet flap of strength Cj and angle e (fig. 2). This 
gives a downwash field w-| . To satisfy the wing boundary condition, a further 
vorticity is required, and this is assumed to be that of a plain airfoil at 
giving downwash Wg . Thus, as shown in reference 1, the kinematic boundary 
condition on the wing is exactly satisfied, and the dynamic boundary condition 
on the wake approximately satisfied. The above model is called the alpha 
effective model. 

An alternative approximation is to assume the vorticity is that 
corresponding to a two-dimensional jet-flap airfoil at the true effective angle 
a - a.- but at a reduced jet strength. The new effective blowing coefficient 
Cj e is obtained by maintaining the same momentum lift in the three-dimensional 
wing as in a two-dimensional jet- flap airfoil 

c j ( 0 + a " a =c) = C j e ( e + 01 " a i) (4) 

This is called the effective model. 

J 

Thus, the loss in jet effectiveness due to three-dimensional flow 
(essentially involving the downwash due to the trailing vortex field) may be 
accounted for either by a reduced effective angle of attack or by a reduced 
effective momentum coefficient. 

Combination of Trefftz and Wing Equations 

It has been shown that introducing a as an unknown requires two Trefftz 
and wing equations to close the set. The Trefftz plane solution gives 
unambiguous results for C n and C t . The wing plane gives two alternative 
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expressions for C n and . Thus, in principle, these are four consistent sets 
which could be chosen. For our purposes, the a effective assumption was 
selected for both C n and C^. since this is slightly simpler algebraically. 


Coefficients Using a Assumption 

For this model, the lift coefficient due to the jet- flap airfoil at ,0 

is C 0 a + C a 0 (where C. and C„ are the lift slopes for the jet- flap 

a 0 a 0 

airfoil at the appropriate C. ). Expressions for C„ and C„ are given in 

J a 0 

reference 2. In addition, there is the lift due to the plain wing vorticity 
of a flat plate at , 2 tt • Thus, the lift coefficient is 

C = 2tt a, + 0 (5) 

a 0 


Subtracting the jet momentum C. (a, + 0) in order to obtain the normal pressure 

J * 

lift coefficient gives 

C t p ’ 2 " “2 + (\ - C j) “1 + (\ - C j) 0 (6) 

For the nose thrust, the singular behavior near the leading edge is required. 

This may be represented in the form of a nose source N as described in 

0 1/2 

reference 5 such that the perturbation near the nose behaves like N Q /x 0 
where x 0 is the distance from the leading edge. It is shown in reference 5 
that this can be written for a two-dimensional jet flap as 


N = N a 
o o 

a 


+ N 0 


(7) 


For the plain airfoil vorticity, N = 1.0. Thus, the total nose 

°a 

singularity constant is given by 


N = N a, + a 0 + N 0 
0 °o 1 2 °8 


( 8 ) 
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As shown in reference 5, this singularity accounts entirely for the nose 
thrust which is given by 

eV (9) 


Jet-flap theory (ref. 5) gives 


C, = 27 r 






( 10 ) 


N 2 = ^ 
°0 4tt 


( 11 ) 


These expressions may now be coupled with the results from the outer solutions. 

Coupling 

Combining the Trefftz and Wing solutions (eqs. (2) and (6) ) for the pressure 
lift coefficients gives 


- Cj (5 + y + o. + 


e ) " 2it a 2 + ( c ^ a " c j) a i + (S 


e J 


c.) e 


Substituting the geometrical results 


ou = a - a = a + a 6. 

1 OO r V 


and 


a-, = o a 
I 00 


( 12 ) 

(13) 

(13a) 


results in 


4 a 


2ir 


(o - 1^ a <f» 


y + C l a (° + 5 *y) 


+ c. 0 

*0 


(14) 


15 



Applying the same procedure to the thrust equations (eqs. (3) and (9))" 

gives 


2lT '{ N o a ( a + “ *y) + N 6 + a <t> y (! - a) } 2 = 



Relationship With Prandtl Equation 

It should be noted that equation (12) is actually an equation connecting 
<f> and <|>y and represents the boundary conditions on the harmonic function 4> . 
These boundary conditions are known as Riemann-Hilbert-Poincare conditions; 
where one part of a harmonic function is defined as a function of derivatives 
of its conjugate. In order to transform this into the classical integro- 
differential form, it is only necessary to note that 




(16) 


by the harmonic property of <}> . Noting that 4> is proportional to the 
circulation r , this is recognized as the classical downwash integral of 
lifting-line theory. Thus, the equation representing the boundary relationship 
between <|> and <f> coupled with the harmonic property of <J> will always give 
rise to an integro-differential equation relating <J> and / (<|> v /t - x) dt . 

SOLUTION TECHNIQUE 


General 

The method employed is to assume a a distribution and use collocation 
methods to solve equation (12). Then the results of this solution are 
substituted into equation (15), and a new result for o is obtained. 
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Normalization of Lift Equation 


Equation (12) can be rewritten in the form 


4 a <f> 



- 2tt (1 - a) 


a + 0 

c &a - 2 tt (1 - a) | + a 


( 17 ) 


Putting Y = 2y/b , X = 2x/b , c* = c/c av , c av = S w /b , A = b/c ay , and 
$ = ci^/c,,, where b is the wing span , S.. the area , and c the mean 

3 V W a V 

chord gives 


4 $ 



A a a + ^0 9 2 

C„ - 2tt (1 - a) l A *Y 


(T8) 


This has the same form as the Prandtl equation for a plain wing of lift slope 2 tt 
euii valent chord c g , and twist a g , where 



(19) 


% a + % 6 



so that equation (18) becomes 


( 20 ) 


4$ 

TT C, 


= 2 


. 2 . 

“e + A $ Y 


( 21 ) 


Potential Functions 

Continuous potential . - Writing the spanwise coordinate X = cosco, the 
continuous potential $ c can be defined as 
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r-1 

^ ■ 12 A n si " n 
n=l 


0) 


( 22 ) 


r-1 

4 - y 

in a) / ^ 


sin 


n=l 


n A n sin n u) 


(23) 


It can be shown for r even (ref. 10) if $ c is defined at the pivotal points 
to = m ir/r with 


r-1 


C = 


m sin 


7 A sin n 2UL 
n a) n r 


n=l 


(24) 


that the gradient $ Y can be written at each pivotal point as 


r-1 


(k) -V $ c 3 

n 


kn 


n*l 


where 


(25) 


1 


^ " r sin !LIL n . cos 


1 


1 


^2 sin n tt^ 


( n + k ) 1 " 1 roc ~ k ) 1 
r i cos r 

for k + n odd 


for k = n 


(26a) 

(26b) 


= 0 


for k + n even, k f n 


(26c) 
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Thus, if a continuous distribution is considered only and equation (21) satisfied 
at the pivotal points, the following is obtained 



where 


s kn = 1 k = n 

(28a) 

s kn ’ 0 k * " 

(28b) 


It is sometimes convenient to compute the potential gradient at other than the 
pivotal point. It can be shown that this is given by 

r-1 

- $y(w) =5^ B*(w) (29) 

n=l 

where 



( 30 ) 
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Discontinuous potential . - For discontinuous wing geometry, it is necessary 
to have a potential function which is continuous, with the proper behavior at the 

wingtips and at infinity and which has a finite jump in its derivative $ iy at a 

★ 

given station w d^ . The function $. which satisfies this has the property 


< (») - l 




sin 

(<* + “d^ ' 

* 

(cos CO - COS ( 0 ^ 

1 i 

\ l . / 

•- + a)^ ^ s i n (jo 

sin 

(<» ~ “d^ 



2 ) J 


$* (go) = - 

’y 

H ( 

o. " 

d i 

“) 



(31) 


(32) 


where H is the Heaviside function with the property 


h (“d, - “) 

= i 

for 

CO < <0 ( 

h (“d, - “) 

= 0 

for 

CO > C0 ( 


(33) 

(34) 


It can be seen in the analysis that this function is capable of exactly matching 
a discontinuity in either c e or a g at Uj . 

General discontinuous solution . - It can be assumed that discontinuities 

in effective geometric data (c g and a g ) occur at s points, given by 

oj. , i = 1 . . . s . Then the general potential can be written in the form 
a i 

r-1 s 

$ (go) = A n sin n a) + A.. $. (co) (35) 

n=f i=l 

i. u 

where is the magnitude of the i tn discontinuity in downwash. 
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Now substituting into the general normalized lift (eq. (21)) at the pivotal 
points, this gives 


r-i 

E 

nsl 



kn 



+ 


+ 



i=l 





(36) 


where 


c is the equivalent chord c Q 
e k e 


at the kth pivotal point and similarly 


with a Q . 
e k 


This yields r - 1 linear equations for the r - 1 + s unknowns , $ n 
and A.. . The additional s equations are obtained by considering the jump 
condition at each discontinuity. 


Jump Condition 

The jump operator on the function g(x) can be defined in the 
following way 


Lim ( ) 

g(x) = 5-H) | ' 9 ^ d + 6 ^[ 


(37) 


Equation (35) can be written 

s 


£ $ C U) + ^ ^ x i M 10 ) 2 a e c e ' A c e f n ^n^ 

i=l n=l 

s 

• l c e E X 1 H ("d, - “)' 


(38) 
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The jump operator is now applied at each discontinuity. Then, noting that 


$ c and are continuous throughout, this gives 


0 = 2 a c I 
e e 'a 


' d .-A c e l» d- £ *n (“d.) 

J j n=l J/ 

s 

' f { c e 2 x i H ( a) d i ' “dj) } 


. . j = 1 . . s 




Now, defining 


the set obtained is 
r-1 


E C (“dj *n + E J f j A i ■ fl A (“e c e)j 

n=l i = l 

j = 1 •’ s 


(3?) 


A(a c ) . = (a. c o ) 1 
e e'j e e' 'u. 

j 

= 1 " s 

(40) 

d j 




- “dj) } l Bdl 

j 

i 

= 1 " s 
= 1 •• s 

(41) 

C (“dj) ■ K ("dj) C e } l* 

j 

n 

= 1 •• s 
= 1 •* r - 1 

(42) 


(43) 


This provides the remaining s linear equations. Equations (36) and (43) 
are then inverted for the solutions and . From these, the potential may 
be evaluated by substituting into equation (35). 
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Auxiliary Thrust Equation 

Once the potential has been determined for the first assumed a distri- 
bution, the thrust equation is checked, and a multidimensional Newton's method 
used to establish the variation of the thrust equation with a . 

The thrust coefficient at the wing plane is given in normalized form from 
equation (15) as 


c ? - 2-J 


(° + ■-$ 


+ N 0 - 2 $ Y 

0 ' 


H) 


(44) 


In the Trefftz plane, the thrust coefficient can be given by the right-hand side 
of equation (15) as 


.t _♦*(•♦$ tc l ° 2 - (a + M 


(45) 


At each pivotal point or singular point, the error function F., is 

K 


F k = c t {k) - c t (k) 


k = 1 ••• r - 1 + s 


(46) 


Then, a at each point can be varied by a fixed value e , generating a matrix 


E 

kj 

The change in a 



h 1 £ ) : F k 


Aa. required to make 

J 



k = 1 •** r - 1 + s 
j = 1 • • • r - 1 + s 

F. zero is now given by 


(47) 


(48) 


23 



NUMERICS AND ITERATION PROCESS 
General 

The procedure used is to solve equations (36) and (43) for 0 = 0.5 , then 
to substitute the values of $ into the auxiliary thrust equation (48) to 
determine a new set of a values to insert in equations (36) and (43) and to 
continue this process until the maximum change in a is less than an arbitrary 
assigned value. 

The program was written with the capacity to handle 8 points of discontinuity 
and 21 pivotal points. The only numerical constants required are the blown-lift 
coefficient slopes. This analysis uses the results given in reference 1: 

Cjt = 2tt (l + 0.151 C. 1/2 + 0.291 C.) (49) 

^ \ J J/ 

C Aq = 2 (tt Cj ) 1/2 (l + 0.151 Cj 1 ^ 2 + 0.139 C^ 1/2 (50) 

It was found that only four iterations were required to bring the interative 
changes a down to 10 - 5 so that convergence in all cases was very rapid. 

Programed on the CDC 6600, a given wing planform takes less than 2 seconds to run. 
The program is written with attention to divisors, so that the case C. = 0 can 

J 

be treated directly. Thus, plain unblown discontinuous wings may also be handled 
by the program. 

Local Integral Properties 

Three items of interest are the local lift, induced angle at infinity, and 
induced drag. If $ is the solution, it can readily be shown that 


= ii. zc *-l 

c* a j A 

(51) 

a = - 2 -f 
00 - A 

(52) 

a 

c d-j = c n 

(53) 
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These coefficients can be integrated across the span (taking due account of 
singularities) to give the integral wing properties and . 

THICKNESS EFFECT 

The preceding analysis is done for airfoils of zero thickness. The thickness 
effect can approximately be incorporated by using a modified lift slope in the 
wing-plane equation, as is normally done in the lifting-line approach to wing theory. 

On this basis, the lift slope of a plain wing of finite thickness can be 

written as 2 it n and that of the jet-flap wing of thickness t as C„ and 

t a 

C„ . Expressions for the latter quantities are derived in Appendix B. Substi- 

tuting these modified lift slopes into the wing-plane lift equation, it can be 
found that the equivalent chord and angle of attack (eqs. (19) and (20)) become 



With these new definitions, equation (21) is unchanged. 

EFFECTIVE ASPECT RATIO 

Aspect ratio plays a dominant role as the scaling parameter in the lifting- 
line equation since the solution depends in a continuous, but nonlinear, fashion 
upon aspect ratio. 

For blown wings, the significance of the geometrical aspect ratio is obscured 
by the additional parameters Cj and 0 . However, it can be shown that, for an 
elliptically loaded blown wing, an equivalent aspect ratio can be defined which 
has many useful properties. 
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The basic integral equation "for the jet- flap wing is 


C« a + C. 0 C - 2 tt (1 - a) 
4 « _ o a *6 . 2 

2 — + T z 


TT C 7 


2tt 


( 56 ) 


Differentiation with respect to a gives 


4 $ 
a 

tt c* 


.-*♦ 2 {C, . 2n (, - 0 )tV 

TT A TT l a ' J 1 


(57) 


An unblown wing, where a = 1/2 and = 2 tt , gives 

a 


If 


4 $ 9 

7T c* 4 A 




(58) 


(59) 


and 


equation (57) becomes 


2 { f -V°)} 


4 ft. 


a 


2 ft 


TT f C* 


= 2 + 


ay 


f A. 


(60) 


(61) 


This converts equations (57) and (58) into identical equations. If f and a 
are not functions of spanwise position, equation (61) can be solved by using 

standard theory for an unblown wing of aspect ratio A g . 
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For the el 1 i pti cal ly loaded case, this assumption of spanwise invariance is 
valid and f can be determined from two-dimensional jet-flap theory. However, a 
cannot be determined without an auxiliary equation. Assuming a is known, where 
is the mean wing-lift coefficient, equation (61) can be solved to obtain 


giving 



2 TT 


n + 2 c j1 

2 

l 1+ /rJ 


i'M 






(62) 


(63) 


where is the two-dimensional lift slope. A similar result can be shown for 

Cl 0 . Thus, the interesting general result that can be obtained is 



(64) 


This result is true for all el 1 i pti cal ly loaded high-aspect-ratio wings - blown 
or unblown. The induced drag may also be written (exactly for an el 1 i pti cal ly 
loaded wing) as 


C Dt = 


I TT A + C- 


(65) 


For a wing of arbitrary loading, f and g will in general vary across the 
span; however, these equations will be approximately valid for some mean values 
of f and a and can be written as 
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where 



1 


1 + f 

H e 


1 + 


2 C J 
ir A„ 



*1 



TT A + 2 Cjj 


A e =A[2 {'f- (1 -.5)}] 


( 66 ) 


(67) 


( 68 ) 


and e is an induced-drag efficiency which will be less than unity. 

It should be noted in this representation a must be known and cannot, in 
fact, be determined without solving the lift and thrust equations. However, 
these results can be useful if an estimate for a is known from previous experience. 

RESULTS AND COMPARISON WITH THEORY AND EXPERIMENT 


Check of Analysis 

It should be noted that the various known limit cases are corre tly obtained 
with this model. Putting the aspect ratio equal to infinity recovers the proper 
two-dimensional jet-flap airfoil performance. For elliptical, loading, a constant 
a across the span is obtained with uniform downwash and lift coefficients. For 
zero blowing, the proper lifting-line solution for the wing is obtained. The 
discontinuous input was found to check with results given in reference 10 for an 
unblown wing with twist discontinuities. 

A severe test of the discontinuity analysis was made by taking an inverse 
problem as a test case. Here, a uniform downwash was assumed, implying an 
ellipticaly loading. However, the wing pi an form was given a chordwise discontinuity 
with the elliptical loading maintained by appropriate discontinuities in angle of 
attack. This particular discontinuous geometry was used as input into the direct 
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program which then inverted to the integral equation to recover the uniform 
downwash and elliptical loading to the fourth significant place. 

A check of linearity was made by determining the lift slope for a set of 
increasing angles of attack from 1° to 30°. It was found that the slope remained 
the same to the fourth significant place. This is interesting in view of the 
fact that, although the lift eqif&tion is properly linear, the auxiliary equation 
for a is actually quadratic in a . 

Basic Results 

There is such a large set of results available from the program that it is 
of no particular value to display wing performance for some arbitrary set of 
geometries. The elliptical cases have already been given in reference 1. 

Thus, for design purposes, the next most generally significant wing could be 
considered to be the constant-chord straight wing with constant spanwise blowing. 

A series of results for varying aspect- ratio and blowing coefficients are shown 
in figure 3 for rectangular wings of 0-percent thickness. 

It is convenient here to plot only the pressure-lift coefficient; that is, 
the wing-system lift after the jet-momentum contribution has been subtracted. 
Typically, at a station, the pressure-lift coefficient is given by 

C SL = C £ - C j ( 0 + «) (69) 

while the wing-pressure-lift coefficient is given by 

C, = 1 f C c db (70) 

P 5 J *P 

It should be noted that c L p is really the significant term in any evaluation 
of the solution technique, since the momentum lift may be mechanically and exactly 
computed for any planform. By removing the momentum lift, any differences in the 
solution approach may be more clearly identified. Because of the linearity of the 
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analysis, the results are further condensed by considering only the gradients 
with a and 0 - that is, and ^Lp 0 . It can be seen from figure 3 that 

this results in a fairly compact presentation of performance. It is also note- 
worthy, that for blowing coefficients above 4, there is little variation in the 
lift slopes with increasing Cj . 

The thickness effect is very significant. For an unblown wing of constant 
thickness, it csm be shown that the thickness term can be accounted for by writing 
2 n n as the section lift slope; then, the lift slope of the thick wing becomes 
n multiplied by the lift slope of a O-percent-thick wing of aspect ratio A/n . 

The effect is to increase the three-dimensional lift slope by an amount somewhat 
less than n . A similar effect occurs for the blown wing; however, this appears 
to be rather large when referred to pressure lift. Figure 4 shows a comparison 
of pressure-lift slopes for O-percent-thick and 1 O-percent-thick wings. It can be 
seen that the thickness is as significant as aspect ratio in determining performance. 
This point becomes important in comparing the theory with analytical results. 

The generalized induced drag efficiency e (eq. (67)) is a useful quantity, 
since this indicates how well a given wing approaches its minimum blown induced 
drag. For rectangular wings, these values of e are shown in figure 5. It can 
be seen that the efficiency values are quite high and suggest that for uniform 
blowing the induced drag of rectangular wings can be approximated by the well- 
known elliptical blown induced-drag expression. 

Comparison With Other Theories 

It should be noted first that since the analysis is an approximation to the 
actual lifting-surface model, which is itself an idealized model of the actual 
flow, the theoretical comparison of existing approaches should initially be 
discussed. The present paper is a lifting-line analysis and, in this respect, 
is similar to those by Maskell and Spence (ref. 1), Kerney (ref. 2), Tokuda 
(ref. 8), and the section design method of Lopez and Shen (ref. 9). The latter 
three methods do not use an auxiliary equation to determine a , but assume it 
has its unblown value of one half. 
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For elliptical loadings, the present model correlates with that of Maskell 
and Spence (ref. 1), since it is actually a generalization of their technique, 
although a slightly different form of the auxiliary equation is used. The proper 
test of the lifting-line/auxiliary-equation model is the comparison with a 
properly numerically modeled linearized lifting-surface theory, with the appro- 
priate kinematic and dynamic boundary conditions applied to a wing of zero thick- 
ness. The only existing lifting-surface model properly satisfying the wake 
boundary conditions is that of Lopez and Shen (ref. 9). Reference 9 shows 
excellent correlation of Maskell and Spence's model with the Lopez and Shen 
lifting-surface approach. The present model can thus be considered verified for 
elliptical loading. 

The rectangular wing presents a more difficult problem of correlation. 
Published solutions are given by Tokuda (ref. 8) and Lopez and Shen (ref. 9); 
these solutions differ by small but significant amounts. As a representation of 
the differences between the various models, table I shows the ratio of the inte- 
grated wing lift to the section lift (C^/C^) for a representative design case - 

elliptical and rectangular wings of aspect ratio 6.0 at the blowing coefficients 
of 1.0 and 2.0. In obtaining these data, the data from reference 9 was ready 
from a graph; while for Tokuda's result (ref. 8), his equation 6.6 was used. 

This equation is read as 


'L . 


1 


1 + 



+ 1.577 


a 

4irA 


(71 ) 


This form is written out since there is some ambiguity in the enunciation of 
equation 6.6 in reference 8. As a further baseline, the classical lifting-line 
unblown (Cj = 0) results of Glauert are shown in the last column. It can be 
seen that the present paper correlates rather well with the lifting-surface 
solution while Tokuda's result (ref.. 8) is somewhat higher. It can be concluded 
that the present analysis correlates with the lifting-surface approach at least 
in this range. 


31 



Table I also illustrates that there is about as much variation in solutions 
between the rectangular and the elliptical planform as there is in the different 
analytical approaches. While all solutions are probably acceptable for preliminary 
design purposes, it would be very valuable to have a few known exact solutions 
to test approximate methods against. 

TABLE I.- RATIO OF WING LIFT TO SECTION LIFT FOR 
O-PERCENT-THICK ASPECT-RATIO- 6 WING 




C T = 1.0 

C T = 2.0 

Cj = 0.0 



J 

J 

(from Glauert) 


Maskell & 

0.699 

0.680 



Spence 
(0 = 0.5) 




Elliptical 

Maskell § 

.700 

.702 

0.751 

Planform 

Spence 
(correct a ) 





Kemey 

.740 

.743 



Tokuda 

0.790 

0.765 



Lopez $ 

.683 

Not avail. 


Rectangular 

Shen 




Planform 

Lissaman 
(0 = 0.5) 

.664 

.656 

0.720 


Lissaman 
(correct a ) 

.675 

.668 



The effect of thickness on lift slope is very significant and appears to be 

larger for the blown wing than a plain wing. Some representative results are 

shown in table II where the thickness correction of reference 3 is used for the 

blown elliptical wing, and the present method used for the rectangular wing. In 

r t r o 

table II, the ratio of u l d for a 10-percent- thick aspect-ratio-6 wing to U L D 

for a O-percent-thick aspect-ratio-6 wing is presented. 
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TABLE II.- EFFECT OF THICKNESS ON PRESSURE-LIFT-CURVE SLOPE 


Planform 


Cj = 0 

Cj = 2.0 

Elliptical 

Rectangular 

1.075 

1.046 

1.100 

1.102 


It can be seen from figure 4 that a 10-percent- thick aspect- ratio-6 wing has 
about the same performance as a 0-percent- thick aspect- ratio-8 wing. 

Comparison With Experiment 

From the preceding comments, it can be observed that comparison with experi- 
ment is a difficult process, largely because few carefully controlled experiments 
exist. Difficulties in proper normalization of experimental data involve not only 
the usual problems of viscous boundary-layer effects which occur with all wings, 
but the special jet- flap problems relating to proper determination of jet coeffi- 
cient and jet angle, as well as the problem of proper wind-tunnel corrections for 
a jet-flap wing. In this light, the approximate thickness corrections and the 
use of an experimental effective aspect ratio can be used to fit a theoretical 
solution to almost any experimental data. 

As an illustration of this point, figure 6 shows a comparison of the present 
theory with the experimental results of Williams and Alexander reported in 
reference 11. The curves of the present theory were obtained by directly inserting 
the given data (namely planform, aspect ratio, thickness, angle of attack, blowing 
coefficient, and jet angle) into the program with no "corrections". It can be seen 
that the correlation is very good for lift slope, but that there appears to be a 
systematic error in the lift at 0° angle of attack. This could be due to the fact 
that the jet is not actually issuing at its nominal angle of 31.3°. It is of 
interest to note that after making an approximate global thickness correction, 
the results of Tokuda (ref. 8) match the experiment even better. However, Maskell 
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and Spence (ref. 1) match the same data with a similar global thickness correction 
although their solution is for an elliptical wing with elliptical blowing, giving 

quite different results for C, /C 0 as shown in table I. 

a a 

CONCLUSIONS 

An analysis has been developed for high-aspect- ratio jet-flap wings of 
arbitrary (including discontinuous) geometry and blowing. This analysis has 
been programed, and the results lead to the following conclusions: 

1. The technique is as accurate as other known theoretical methods. 

2. The program is more flexible and rapid than other methods, 
particularly with regard to discontinuities in wing geometry. 

3. Correlation with existing experimental data is satisfactory, but 
it is considered that definitive test data are lacking. 

4. The present method is a useful analytical tool for jet-flap wing 

design. 

5. Two additional steps are needed to evaluate jet-flap theories. 

These theories should be checked against: 

a. Definitive lifting-surface analyses 

b. Properly controlled and corrected experimental tests 

Such checks should include cases with variable discontinuous planform and 
blowing parameters. 


34 



APPENDIX A 


TREFFTZ PLANE FORCES 
GENERAL 

Here, the situation in the Trefftz plane where all forces on the wing may 
be indirectly determined is considered. It is convenient to use the exact 
nonlinear large angle analysis initially, and then to expand to the appropriate 
order of accuracy. 


UNBLOWN WING 

The unblown case is treated first to establish a reliable connection with 
conventional Trefftz plane theory. It is assumed that a self-preserving steady 
wake exists in the Trefftz plane (fig. Al). Thus, in places inclined at a to 
the mainstream flow V , the situation is properly two dimensional. Incompressible 
flow of density p is assumed. In this inclined plane, for steady two-dimensional 
flow, a complex potential W(z) is defined as 

W = V sin a (- iz + w(z)) (Al) 

where z = x + iy with x horizontal and y upwards in the Trefftz plane. 

The normalized perturbation potential w(z) may be written 

w(z) = <j> (x,y) + i ¥ (x,y) (A2) 

Forces in the direction z and y are now considered by applying the 
momentum theorem to a large cylinder of axis parallel to the final wake direction. 

On the upstream face, velocity and pressure have their free-stream values; 

on the downstream face, the velocity is V cos a in the direction of the 

cylinder axis V sin a <t> Y in the x (spanwise) direction and V sin a (1 + <f> ) 

x y 
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in the orthogonal y direction. Thus the pressure perturbation Ap in the 
Trefftz plane is given by 

Ap = j p V 2 1 1 - cos 2 a - sin 2 a (1 + «J> - sin^ a <J> X 2 j (A3) 

For the forces in the Z direction F z , there are only pressure 
contributions giving 

F z = -f Ap dS* (A4) 

S 1 

where S-| is the infinite downstream face of the cylinder and 
elemental surface area, substituting equation (A3) gives 

2 F z = p V 2 sin 2 of ^2 <J>^ + <j> x 2 J dS* 

S 1 

For the y force, only momentum need be considered giving 

F = -p V 2 cos a sin a f cf> dS* (A6) 

y j y 

s i 

Noting that 4> decays sufficiently rapidly at x,y ». , these surface 
integrals can be converted by Gauss Theorem into line integrals about the wake. 
Assuming the wake is linear and oriented in the x direction, the symmetries of 
the upper and lower surfaces of the wake give 

2 F z = p V 2 sin 2 aj-4y~<j>dx-2 <J>^ dx j (A7) 

Fy = p V 2 cos a sin a 2 .J'fy dx (A8) 

where <j> is the potential and <t> y its gradient on the upper surface of the wake. 


dS* i s the 

(A5) 
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This can now be resolved into the conventional lift and drag directions to 
obtain the following exact relations 

L = p sin a | 2^$ dx + sin^ afty <j>y dx} (A9) 

D = - p sin^ a cos a^J) <J> y dx (A10 ) 


For a small, these degenerate to the classical linearized results. Results 
similar to the above are given by Ribner (ref. 12); however, in that version, the 
induced drag is left in the form of an infinite surface integral, as opposed to 
the convenient finite line integral given here. 


It is now assumed that the trailing vorticity which must satisfy the Beltrami 
condition, that is, be parallel to the local flow, has negligible lateral 
(spanwise) components in the wake. Then, each spanwise station in the Trefftz 
plane contains the same lift as the corresponding spanwise wing station. Thus, 
defining the local lift coefficient as C^(x) and local chord as c , the 
following is obtained 


C^(x) = 4 sin a 


4(x) + sir> 2 ■ P - <fr(x) <fr y (x)} 


(All) 


By making the same assumption for the local induced-drag coefficient (x ) , the 
following is obtained 


2 - - T v 

C ( j(x) = - 2 sin a cos a 4>(x) - JL 


<k,(x) 


(A12) 


It should be noted here, that while this result is correct for the lift for 
any planform, it is only valid for induced drag if the planform is essentially 
unswept. For swept wings, it is known that the distribution of induced drag in 
the Trefftz plane does not correspond to that on the wing itself, although the 
total induced drag, of course, can be obtained by integration either on the wing 
or in the Trefftz plane. 
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BLOWN-C-ASE- 


The jet sheet is now introduced assuming the momentum at each station to be 
y(x) and inclined at an angle 3(x) to the horizontal (positive downwards). Then 


3 = a - tan~^ j tan a (1 + 4> y )j 


3 a - a <j>y + a (a 3 ) 


(A13) 

(A14) 


Thus, to second order, the additional force coefficients due to the jet 
element are 




2 y(x) a ♦ (x) 
“c^ 


(A15) 


c d j(x) = - 2 Y (x) 


(l - | ^(x) 2 ) 


c p 


V2 


(A16) 


Then, writing C- as the local jet-momentum coefficient with c the local chord 

J 




(A17) 


for the blown wing gives 


c ^(x) = + 4 ”^( x ) - Cj a <D y (x) + 0 (a 3 ) 


(A18) 


2 a 2 <|>(x) «fr (x) ( a 2 d> v (x) 2 \ /,v 

C d< x > c * - C j l 1 fc— j + 0 (« > 


(A19) 
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It is shown in equation (A19) that the terms without C- account for the vortex 

J 

lift both on the wing and in the jet wake. 

The force contribution due to the entire jet sheet can now be computed 
directly by considering conditions at the trailing edge of the wing. Assuming 
the wing to be uncambered, with the local geometrical angle of attack to be a(x) 
and with the jet angle relative to the wing to be 0(x) , the direct momentum 
force is subtracted to obtain the pressure components on the wing itself in 
coefficient form as 


C ^p = + 4 "c ^ " C j (“ ^y + a + S ) ( A2 °) 

S = ^ c j {a 2 4> y 2 - (<* + e) 2 } ( A21 ) 

These are now resolved relative to the chordline to give the normal force 
coefficient C n and the thrust coefficient C t . To second order in a and a , 
these become 


C 


n 



+ a + 


6 ) 


(A22) 


C t = ^4^ ( a + + j® 2 ( a + 5 ^y) 2 } (A23) 

This pair constitutes the outer solution, expressing the pressure forces on the 
wing itself. 
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JET-FLAP THICKNESS CORRECTION 
GENERAL 

The analytical results for jet-flap performance used in the paper are for 
the linearized problem; that is: the airfoil thickness is linearized; the position 

of the jet wake is assumed to be on the airfoil chordline; and the dynamic boundary 
condition itself is linearized. To find the first order effects of linearization, 
an exact nonlinear solution is assumed to be known for a Joukowski airfoil, and 
they by a mapping to an auxiliary plane. This is reduced to a flat-plate airfoil. 
Then, by expanding the mapping to appropriate order, the ratio of the lift in the 
two planes can be determined, assuming that the lift on the flat airfoil in the 
auxiliary plane is given with sufficient accuracy by linearized theory. The 
analysis is two dimensional and incompressible. 

ANALYSIS 


Figure B1 shows a thick Joukowski airfoil and the auxiliary airfoil and 
defines the symbols used. The potential in the physical plane, Z = X + i Y given 
by W(Z) = $ + i V , is assumed. Q is the complex conjugate velocity dW/dz . 

The boundary conditions are as follows: 


kinematic: *F = 0 on airfoil and wake 

2 

I / A\ 9 / 9 I ^ ^ ^ 

dynamii 


*c: [(Q + ) 2 -(Q-) 2 ]=^ 


(Bl) 

(B2) 


where + and - refer to upper and lower surfaces of the vortex wake , C the 
chord, C- the momentum coefficient , and R the radius of curvature of the wake. 

This geometry is now assumed to be mapped to the auxiliary z' ° x + i y 
plane by 
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z' = F(Z) , dz'/dZ = f* (B3) 

where under f* the airfoil maps to a slit. Potential is conserved at 
corresponding points. The mapping is required to have the following properties: 

f* + 1 , Z -*■ « (B4) 

f* regular at T 1 (B5) 

fi? analytical elsewhere except at N (B6) 

where T 1 is the trailing edge of the airfoil and N the leading edge. It should 
be noted that corresponding elements are given by 

dz' = f * dZ (B7) 

so that the length and angle of an element in each plane are given by 


ds = |f*| dS , 9* = (arg f*) + 0* 


(B8) 


with 


dz 1 = ds 6 10 * , dZ = dZ <J> 19 * 


(B9) 


Under this mapping, the boundary conditions become 


kinematic: ¥ = 0 on body and wake 


dynamic: (q + ) 2 - (q‘) 2 J = c‘ q ro 2 ~t -ft 


I 1 _ d!* 7 (arg f *)} 3s“ •• on wake 


(BIO) 
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where 6* is the inclination of the wake in the auxiliary plane; ds an arc 
length along the wake in the plane; and q the complex conjugate velocity in 
this plane. 

Thus, this now becomes a nonlinear problem of flow about a jet flapped flat 
plate of chord c 1 at angle a with jet angle of 9 and with a momentum 
coefficient c^ which varied according to 

C J - 35* f *>} < B ”> 

THE MAPPING FUNCTION 

The parametric functional relationships can be expressed as 

Z = y + 1 (B12) 

r 


2 ' 


y + a + 


Ltii 


(B13) 


where y is an intermediate complex parameter. This transforms the Joukowski 
airfoil to the flat plate, with the required properties on the derivative. In 
addition, the following geometric properties are noted 


Z Plane 


Chord = 4 + — a ' 


1 + 2 a 


= C 


z' Plane 


4 (1 + a) = c 1 


Thickness 2 y - + - - a 


= T 


0 = t 


The mapping derivative 





dZ 


(B14) 
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is evaluated as 


f* 


I * “ n 


► < 

*• 

— M 2 . 

J 

y 2 - 1 


^1 - 2a 


/ 1 ~l " T + 0 ( a3 ) 
y(y + l) v 


At the trailing edge 


f* <v 1 - a - a^ 


while downstream 


f* 1 


(B15) 

(B16) 


(B17) 


(B18) 


The approach to unity is very rapid. For example at y = 3 , which is about 
0.5-chord length downstream, it can be found that f* = 0.98 for a 10-percent- 
thick airfoil . 


EQUIVALENT PROBLEM 

Thus, the equivalent problem consists of finding the exact solution in the 
z plane of the flow about a flat plate of chord c 1 =4 (1 + a) immersed in a 
distance flow of velocity q^ and angle a and having a jet flap of variable 
strength. This strength is given by equation (Bll), substituting equation (B16) 
and expanding to 0(a) gives 


c 


j 



at the trailing edge 


(B19) 


C j 



at infinity 


(B20) 
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Substituting for 


C 


4 (1 + a il 

(1 + 2 a) 


(B21) 


and expanding to 0(a) gives 



at the trailing edge 


(B22) 


c • = 7 T — “rr- at infinity (B23) 

J V " 7 

The pressure field on the airfoil is now assumed to be a function of the 
blowing coefficient Cj-j at the trailing edge. Then, using standard linearized 
theory, y is the circulation about the airfoil alone and is given by 

a 

2 

-ir ■ ( c l a « + c 2 e e ) - (“ + e ) C J1 < B24 > 

where C° and C° are the lift slopes for 0 thickness computed for a blowing 
coefficient of c.^ . 

Spence (ref. 2) gives 

C *a = 2 f C J') V2 (' + ° 

but 

c ji = c j (’ + a2 + 0 ( a2 )) = c j 
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So 



near the trailing edge, giving 



(a + 9) C. 

J 


(B27 ) 


since 

c ° t = c l° + % s (B28 > 

Under the mapping, Z z , potential is conserved so tha circulation about 
the thick airfoil in the physical Z plane , r_ , is given by 

a 

2 r a = C {c° - (a + e) c.| (B29) 

However, the wake lift in the physical Z plane is given by 

(a + 9) C. C (B30) 

J 

Thus, the total circulation in the physical plane is given by 

2 T = C (C° - (a + 6) Cj C) + (a + e) C (B31) 

So that the lift coefficient in. the physical plane is given by 



(B32) 

noting that 


= 1 + a + 0(a 2 ) 


(B33) 
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APPENDIX B 


Noting that the mapping gives 


T a * 1 + 2a 
C = 1 + a * 1 + a^ 


(B34) 


it can be shown that 



(B35) 


This gives the result to 0(T/C) where T/C is the thickness to chord ratio 



C. (a + 0) ) 

C° \ 


(B36) 


Thus, in physical terms, it is observed that the thickness affects only that 
portion of the lift directly experienced on the airfoil (the pressure lift). 

This result was intuitively proposed as a correction by Spence. The above 
analysis, while still an approximation, provides a rationale for this correction. 


i 
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FIGURE 2. CHORDWISE DOWNWASH OF JET FLAP WING 
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